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SUMMARY 


Values for the stability and control derivatives of a single-engine, low-wing, 
general aviation airplane have been determined from flight data. Lateral and longi- 
tudinal transient maneuvers were analyzed by the equation error and output error 
methods. One quasi-steady maneuver was also investigated. An angle of attack range 
between 4 and 26 degrees was covered and there was a good agreement between the 
parameters extracted from flight data and those predicted by wind tunnel. 


INTRODUCTION 


There have been a number of previous attempts to determine stability and control 
parameters from flight data. Gerlach and Howard in references 1 and 2 have applied 
the equation error method for the determination of longitudinal and lateral 
parameters, respectively. In references 3 and 4, parameters estimated from the 
output error method were compared with derivatives obtained from the wind tunnel and 
theoretical predictions. The parameter estimation problem is reasonably routine in 
flight conditions for which the aerodynamics are linearly related to the response and 
input variables. However, problems persist for the modeling of unsteady and high 
angles of attack aerodynamics. Most of the recent research effort has, therefore, 
been directed primarily towards these unresolved problems (refs. 5 to 11). Even 
though the problem of parameter estimation at low angles of attack has been well 
addressed in the past years, there have not been many applications in high angle of 
attack regimes. Furthermore, few attempts have been made to correlate the parameters 
extracted from flight data to the wind tunnel and theoretical predictions. 

This work was done under the General Aviation Stall/Spin Program at NASA Langley 
Research Center. The primary objectives of this program were to gain a better 
understanding of the aerodynamics in pre- and post-stall regimes and to design and 
validate airframe modifications to improve stall/spin characteristics. Part of the 
overall program was devoted to the measurement of the airplane transient maneuvers 
for the subsequent extraction of the stability and control derivatives. The test 
which produced the data for this analysis was tested with several modifica- 
tions both in flight and in the wind tunnel. For the purpose of this report, data 
from only one of these modifications was used. In flight the airplane was excited 
from steady state flight at different airspeeds by conventional control surfaces, so 
that the entire operating angle of attack range was covered. The mathematical model 

was so postulated as to include nonlinear contributions to the equations of motion at 
high angles of attack* 

The purpose of this report is to estimate the aerodynamic parameters from flight 
data covering angles of attack between 4 and 26 degrees, and to compare these 
parameters with those estimated from the static and oscillatory wind-tunnel tests* 
Also the validity of the parameters is checked by comparing the simulated data which 
is generated from the estimated model with the actual flight records. The estimation 
in this report was done by the stepwise regression method. Experience has shown that 
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the stepwise regression and maximum likelihood methods give similar results. 
Therefore, only a few runs were analyzed by the maximum likelihood method to check 
the regression estimates. 


SYMBOLS 


A 

a„ 


m 


"N 


T 


defined In Appendix B 
longitudinal acceleration, g units 

lateral acceleration, g units 

vertical acceleration, g units 

span, m 

axial-force coefficient, -F /qs 
drag coefficient, D/qs 
lift coefficient, L/qs 
rolling-moment coefficient, Mj^/qSb 
pitching-moment coefficient, M^/qSc 
yawing-moment coefficient, M^/qSb 
normal-force coefficient, -F^/qS 
thrust coefficient, T/qS 
longitudinal-force coefficient, F /qS 
side-force coefficient, Fy/qS 
vertical-force coefficient, F„/qS 

La 


C 

D 

G 

g 

h 

^XZ 

J 

k 

Mx.M^.Mz 


m 

N 

P 

q 

q 


mean aerodynamic chord, m 
drag , N 

force along X, Y, and Z body axis, respectively, N 

airplane response vector 
acceleration due to gravity, m/sec^ 

distance of CG aft of leading edge of wing, percent of c 
moment of inertia about X, Y, and Z body axis, respectively kg-m^ 

product of inertia, kg-m^ 

cost function 

normalized frequency, o)c/2V or o)b/2V 
lift, N 

moment about X, Y, and Z body axis, respectively, N-m 
mass, kg 

number of data points 

body axis roll rate, rad/sec 

body axis pitch rate, rad/ sec 

dynamic pressure, l/2pV 
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weighting matrix 

body axis yaw rate, rad/sec 

wing area, 

thrust, N 

velocity along X, Y, and Z body axis, respectively, m/sec 

airplane total velocity, m/sec 

aerodynamic force or moment coefficient 

angle of attack, rad or deg 

angle of sideslip, rad 

defined in Appendix B 

aileron deflection, rad 

elevator deflection, rad 

rudder deflection, rad 

thrust vector angle with respect to longitudinal body axis, rad 

stability and control derivative vector 

pitch angle, rad 

air density, kg/m^ 

roll angle, rad 

yaw angle, rad 

frequency, rad/sec 


Subscripts : 


o trim value 

M measured value 

WT wind tunnel 

Superscripts : 


T transpose 

• derivative with respect to time 

» derivatives defined in Appendices A and B 

estimate 

Abbreviations : 


CG center of gravity 

ML maximum likelihood 

SR stepwise regression 


Derivatives : 
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Flight Test and Data Reduction 

A single-engine, low-wing, general aviation airplane with a fixed tricycle 
landing gear was used as the test vehicle. Since this airplane had been tested in 
spins, a spin recovery parachute system was installed. Furthermore, the 
configuration flight tested for parameter extraction had as an additional 
modification, an outboard leading edge droop, which enabled trimmed flight at high 
angles of attack. Figure 1 is a three view drawing of the test airplane and figure 2 
is a drawing of the tested configuration. Table I lists the physical characteristics 
of the airplane. The flight data was recorded by an onboard instrumentation system, 
and was sampled at 40 cycles/sec. However, for the estimation purposes every other 
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sample was used. For the longitudinal maneuvers, a simple elevator doublet or a 
combination of several doublets was chosen based upon reference 4. For lateral 
maneuvers the aileron and rudder inputs should have not only good harmonic content 
but also proper phase relationships. In accordance with the conclusions in reference 
4, a rudder followed by aileron input or vice versa was chosen. Typical time 
histories for the transient longitudinal and lateral maneuvers are given in figures 3 
and 4, respectively. In addition, one slow acceleration/deceleration run was also 
analyzed as quasi-steady maneuver, where airplane angle of attack is varied gradually 
by the elevator maintaining a constant pitch-rate. 

The validity of the flight data was checked using the technique presented in 
reference 12, which gives a way of estimating the states via the aircraft kinematic 
equations. In case the integrated accelerations do not match with the flight data, 
appropriate biases are computed. Wherever needed such instrument biases and dynamics 
were taken into account. The measured angles of attack and sideslip were corrected 
for the upwash and sldewash. The measured velocity was corrected based upon the 
calibration derived from flying the test airplane over a measured ground trace. 
Since the equations of motion used in the analysis are referred to the airplane 
center of gravity, the data were transformed accordingly. 

Parameter Extraction Techniques 

There are several methods for the estimation of airplane stability and control 
parameters. Their basic differences are in the optimization criteria and assumptions 
regarding external disturbances and the presence of measurement noise in the data. 
The well established methods for airplane parameter estimation are the equation error 
and maximum likelihood methods. The present report uses a stepwise regression method 
which is a version of the equation error method. 

(a) Stepwise regression (SR); This method is outlined in reference 14. It is 
based on the minimization of the cost function 

J(0) = E - y(0)i]^ (1) 

In the equation, y,^ is the aerodynamic force or moment coefficient based upon the 
measured states and y(0) is calculated from the estimated parameters (0). The 

difference between stepwise regression and ordinary linear regression is in the capa- 
bility of the first technique to select Independent variables from a set of candidate 
parameters one at a time, until the regression is completed. The order of insertion 
and the adequacy of the model is determined by various statistical criteria. For a 
more detailed description the reader is referred to references 13, 14, and 15. 

(b) Maximum likelihood (ML) ; This is a nonlinear estimation technique and 
therefore requires an iterative solution. The method minimizes the error between the 
measured and predicted outputs. For no process noise the simplified cost function 
that has to be minimized is 



where G is the airplane response vector, [V a q e]"^ or [3 p r <|)]'^, and R~^ is the 
weighting matrix. 
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For a detailed discussion the reader is referred to references 13, 16, and 17. 
Results and Discussion 


There were 7 test flights for two CG locations, out of which a total of 62 
longitudinal and 58 lateral maneuvers were available for analysis. In addition a 

total of 6 slow acceleration-deceleration runs for the two configurations were also 
made. The stepwise regression method was considered more suitable than the maximum 
likelihood method for the estimation due to its simplicity and its capability of 

selecting the appropriate model. Also, the low noise level on the data guaranteed 
the efficiency of the technique (ref. 10). 

Typical results for the longitudinal and lateral modes are shown in Table II. 
Figure 5 Illustrates an example of a regression model and the autocorrelation 
coefficient of the residuals of the computed and measured coefficients. The features 
to note are the reasonable fit for the computed and measured coefficients, and the 
rapidly falling autocorrelation of the residuals signifying a random sequence. The 
estimated parameters using two parameter estimation methods are shown in figures 6 
through 13 for idle power conditions. The standard errors were under 10 percent for 
most of the parameters. The model determined was linear until about 10® angle of 
attack, and at higher angles of attack nonlinear terms (Appendix A) were also 

Included. The nonlinear terras improved the fit to the data, but did not 

significantly change the values of the linear terms. Since different nonlinear terms 
were selected as being siginificant for each parameter extraction run, no attempt has 
been made to document them. The figures also show the repeatability of results from 
various flights and their consistency with the ML estimates. Some maneuvers for 
different power settings were also available, but the data were insufficient to make 
any useful observations . 

The acceleration/deceleration run was analyzed assuming quasi-steady flight 
(constant pitch-rate). Using the known elevator effectiveness, the effect of the 
elevator was removed to obtain the aerodynamic coefficients as a function of only 
angle of attack. Figure 14 shows a good agreement between the wind tunnel and flight 

determined and C . 

L m 

Figures 15 through 22 give a comparison of the parameters estimated from flight 
using the SR method and the wind tunnel predictions. 

1. Longitudinal Parameters: 

The X-force derivatives generally had higher standard errors (10 to 50 
percent) and had a poor comparison with the wind-tunnel results, as shown in figure 
15. This is attributed to thrust effects (l.e., idle power is not necessarily zero 
thrust) and the lack of excitation of this mode. A definite trend is apparent until 
about 12® angle of attack, and a scatter is observed for higher angles of attack. 

The Z-force derivatives were well defined (standard error less than 5 percent) 
and the static derivatives showed good agreement with the wind-tunnel estimates 
(Figure 16). , however, did not agree with the wind tunnel, which predicts a 

q 

positive trend in the parameter. A study of the aerodynamic characteristics of 
general aviation airplanes in the same class substantiated the results obtained from 
flight (reference 18). The unusual contradiction between the flight and wind tunnel 
estimates has not been explained. 
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■? derivatives were generally of low standard error (less than 

5 percent) and distinct trends with respect to angle of attack were observed as 

own in Figure 17. Though the configuration flight tested had different CG 

locations, the flight estimates failed to separate the respective s (figure 8). 

An attempt to separate the s was made using the ML method, but there was no 

significant change. Overall, the trends between the flight and wind tunnel C s 

m 

agreed very well, c; was harder to estimate due to the high correlation with C “ 

q m^ » 

ML technique. The values estimated were of a smaller magnitude 

othL Sir .rrf ^he wmd tunnel estlmatT on the 

other hand predicted an unusual decrease in with increasing angle of attack, 

which essentially implies an Increase in short period damping. The coraDarlson 
e ween flight and wind tunnel s was therefore poor. The elevator effectiveness 

q 

(C^^ ) was lower than what wind tunnel predicted. A study of the spin parachute 
e 

dlscrepLc^."*"^” ^ possible explanation gave no clue as to the cause of this 
2 • Lateral Parameters : 

The 3-derlvatlves were identifiable (standard errors less than 5 nercent) 
an except or agreed with the wind-tunnel results (figure 18). C. was consis- 

P ^ g 

tent until about 8° angle of attack and then there was an under-estimation of the 
parameter magnitude from there onwards. This could be due to the flow separation 
over the wing planform, resulting in a sluggish dihedral effect. ^ 

The p-derivatlves with an exception of were consistent with the wind-tunnel 

unS'y P^rcen”" standard errora »ere 

p cent. is a less significant parameter since the airplane motion is 

critical^. therefore, the lack of consistency is not considered 

flfiurr2o'''^^qtrnH^^rr Sen^ally consistent at low angles of attack, as shown in 

Thf errors of the estimates were sometimes as high as 20 percent 

The disagreement at higher angles of attack could be due to the spin recover^ 
parachute canister which protrudes about 1/2 meter aft of the elevator trailing edge! 

verv Imall^n eL easily determinable and errors of estimation were 

very small (less than 5 percent). The aileron effectiveness (C, ) agreed with the 

wind tunnel as shown in figure 21 . ^ 

The effect of the rudder also agreed with the wind-tunnel results, standard 
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errors of the estimates were less than 5 percent. A sharp drop in the magnitude of 
Is apparent around 12" angle of attack (figure 22). This could be due to the 

0 

r 

effect of the wing wake on the vertical fin. At higher angles of attack the vertical 
fin is below the wake resulting in an Improvement in the rudder effectiveness. 

An high angles of attack the SR technique yielded numerous nonlinear terms 
(Appendix A), which have not been discussed. The discrepancies in flight and wind 
tunnel determined parameters in such regimes can be attributed to the effect that has 
been absorbed In such terms. 

The ultimate test for the validity of the model was a comparison between the 
measured and simulated response of an airplane when subjected to the same input. The 
simulated data were generated using the parameters estimated from the SR algorithm. 
Models estimated at different flight conditions were tested and in most cases the 
prediction capability was good. Typical comparisons for the longitudinal and lateral 
maneuvers are given in figures 23 and 24. 


CONCLUDING REMARKS 

A complete set of stability and control parameters for a general aviation 
aircraft was obtained from flight data. The standard errors of the main derivatives 
varied between 1 and 10 percent, errors were higher for parameters of lesser 
importance. Most of these parameters agreed with the wind tunnel estimates. 
Additional work is required to explain the discrepancies observed in the longitudinal 
rotary derivatives. Based on the estimates the airplane response was predicted well 
for the 4 to 20 degrees angle of attack range. 
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appendix a 


Equations of Motion 


The airplane equations of motion are referred to the body axes (see figure 25) 
These equations are based on the following assumptions : 

1 . the airplane Is a rigid body, 

2. the effect of spinning rotors are negligible, 

3. the airplane has a plane of symmetry XZ, 

4. the airplane motion from Initial reference conditions consists 

of small perturbations. 


Based on the above assumptions, the equations of motion are expressed as 


u = -qw + rv - gsln0 + 


2 

PV S 

115“ 


Sc 


V = -ru + pw + gcos0sln<|) + 


PV^S 


w - -pv + qu + gcos0cos<|) + 


2m '^Y 
2 

PV S 


2m '^Z 


^Y-^Z 


XZ 


pV^Sb 


^ ^ (P‘1 + s 


^z-'x 


q = pr 


^XZ pV^Sc 

+ — (r^ - p2) + 


^X ^Y 


Y / '^Y 


2Iy m 


r = pq 


XZ 

.-h ^ 

z / 

0 = qcos4> - rsln(j) 

= p + (qsin<|) + rcos<j>) tan9 


PV^Sb 


+ — (p - qr) + — - — c 

1- 2I„ n 


where 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(A8) 


= C^cose + Cj^slna - C^cosa 


(A9) 


C 


Z 


Cosine + Cj^cosa - C^slna 


(AlO) 


For the equation error method, the aerodynamic 
the measured quantities as follows: 


coefficients were calculated from 


mg 

fs ®X 


(All) 
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C = a 
Y qS Y 

C = -^ a 
qS Z 




X 

qSb 


^Y~^Z 




m qSc 
I, 


q - 


'I -I ' 
Z 


pr - 


C = 
n 


_L_ 

qSb 




pq - 


XZ 

^Y 


z '^2 2 . 

(r - p ) 


(p - qr) 


(A12) 

(A13) 

(A14) 

(A15) 

(A16) 


The aerodynamic coefficients were postulated as functions of the state and input 
variables and their combinations : 


(a) The longitudinal coefficients Cx, Cy, and (L, as functions of a. q. 

aq, a6^, 3^ a3^ oJ , a®. 

(b) The lateral coefficients C^, and C as functions of 3, p, r, 6 , 6 

a r 


a3. up, ar, a6^, a6^, a^p, a^r, g2 ^ g3 ^ ^5^ g3„2^ p3^^ 

2 3 

a , a . 

All the above variables and their combinations are the Increments with respect 
to their trim values. Typical linear models for the longitudinal and lateral 
aerodynamic coefficients are as follows: 


^^o ^Z 2? S. ^'^e"'^e ^ 

o a q 0 o 

e 

C=C +C (3-3) + C ^+C — 

" "o "6 ° "p 2V ^ \ 2V 

+ C (« -6, ) + C (6-6 ) 

3 n. ' r r 

0 O 0 o 

a r 


(A17) 


(A18) 


To avoid Identification problems, the linear pitching moment coefficient is defined 

23 C • 


as: 


C* (a-a ) + C* + C' (6-6 ) 

o ® ®- 2V m. '' e e 


m 2V 

q 


(A19) 


where 


C = C + C 
m m 


+ C c. . S£. cose 


o -o K\ ■ 2v2 


(A20) 
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appendix b 


Wind Tunnel Data 


A static force investigation was conducted on a full-scale model of the test 
airplane (reference 19). Various wing and tall modifications were tested for 
determining their effect on high angles of attack aerodynamics. The investigation 
covered angles of attack ranging from -9 to 41 degrees at a Reynolds number of 
2.5 X 10°, based on the mean aerodynamic chord. 

All longitudinal forces and moments were presented in the wind-axes system, 
whereas all the lateral-directional forces and moments were in the body-axes system. 
Therefore, all the longitudinal forces were transformed into the body axes as 
follows: 


Cjj = -CpCosa + C^sina (gj) 

= -CpSina - C^^cosa (B2) 

The derivatives with respect to the angle of attack were determined by putting a 
second order spline through the data and calculating the slope at each point 
(reference 20). 

Based on (Bl) and (B2), the derivatives in the body axis system were calculated 
from the expressions, 

= (Cj^ - Cp )coso + (Cp + )slna (83) 

^2 ~ ^ C^)cosa + (Cj^ - )sina (84) 


The derivatives with respect to the control deflection were determined based on the 
coefficient changes over the entire angle of attack range. 

(B5) 

where A is C^, C^, C^, 0^^^, or C^; and 6 is 6^, 6^, or 6^.. 

Transformation equations for the control derivatives for and follow from 
(Bl) and (B2) as follows 




COSO + sina 

(B6) 

e 

e 

e 




sina + C cos a) 

(B7) 

e 

e 

e 



C was transformed to the airplane Center of Gravity by 
a 
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The values for , C , and C were read directly from reference 19. 

3 3 a 

The dynamic derivatives were derived from the unpublished wind-tunnel tests that 
were run on a one-third scale model. A power spectral density analysis for both 

longitudinal and lateral transient maneuvers was required to determine the dominating 
frequencies. ® 

For the longitudinal case k was determined to range between 0.028 and 0.04 for 
angles of attack ranging between 4 and 24 degrees. A similar analysis on lateral 
transient maneuvers between angles of attack 6 and 26 degrees yielded a k value 
ranging between 0.15 and 0.24, the wind tunnel data however unlike the longitudinal 
case was Insensitive to this variation. 

The following derivatives given In an unpublished wind-tunnel test report are; 


It 

Cy + Gy, slna 

P 3 

(B9) 

c = 

n 

P 

C + C slna 

"p ^3 

(BIO) 

r' =» 
P 

S + C slna 
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C + C 
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q “ 

(B12) 
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q a 

(B13) 
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q « 
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CL - CL, cos a 
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TABLE I.- PHYSICAL CHARACTERISTICS OF TEST AIRPLANE 


Wing (Modified NACA 642-415) 

Span, m 

2 

Area, m 

Mean aerodynamic chord, m 

Aspect ratio 

Dihedral, deg ...... 

Wing incidence, deg . . . 

Aileron (each) 

Span, ra 

2 

Area, m . 

Chord, m 

Flap (each) 

Span, m 

2 

Area, m ....... 

Chord, ra 

Horizontal tall (NACA 65^-012) 

Span , 

2 

Area, m 

Mean aerodynamic chord, m 
Incidence, deg 

Elevator 

Span, m 

2 

Area, m ....... 

Root chord, m .... 

Tip chord, m 

Vertical tall (NACA 65^-012) 

Span, m . . , 

2 

Area, m . . . . 

Root chord, m 

Tip chord, m ....... 

Rudder 

Span, m 

2 

Area, m ........ 

Root chord, m 

Tip chord, m 


7.45 

9.21 

1.23 

6.10 

5.00 

3.50 


1.16 

0.24 

0.21 


1.15 

0.25 

0.21 


2.34 

1.96 

0.84 

-3.00 


2.34 

0.77 

0.34 

0.21 


1.25 

1.20 

1.09 

0.51 


1.25 

0.33 

0.34 

0.21 


Propeller diameter, m 
Propeller pitch, m . , 


1.80 


1.17 
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TABLE II.- ESTIMATE AND STANDARD ERRORS OF PARAMETERS 


Longitudinal 


Parameter 

Value 

Standard Error 


-0.0806 


P c 

0.396 

0.0076 

q 

1.19 

0.263 

e 

in 

• 

0 

1 

0.0109 

V 

1.74 

0.116 


Parameter 

Value 

Standard Error 

S 

o 

-0.0554 


S 

a 

-3.43 

0.0117 

S 

q 

-6.07 

0.406 

% 

e 

-0.214 

0.0169 

2 

a 

2.4 

0.178 


Parameter 

Value 

Standard Error 

c* 

m 

o 

0.00279 


C’ 

m 

a 

-0.223 

0 .00367 

C’ 

m 

q 

-6.21 

0.127 

C’ 

-0.607 

0.00526 

"6 


e 




TABLE II.- ESTIMATE AND STANDARD ERRORS OF PARAMETERS (concl’d) 

Lateral 


Parameter 

Value 

Standard Error 


0.024 


o 



S 

-0.328 

0.0042 


0.519 

0.0204 

p 




-0.374 

0.0505 

r 



a 

0.181 

0.00442 

r 

0.0593 

0.0039 


Parameter 

Value 

Standard Error 


-0.00113 


o 




-0.0827 

0.00174 

c. 

-0.382 

0.00844 

P 




0.214 

0.0209 

r 



S 

6 

a 

-0.108 

0.00183 

r 

0.0152 

0.00162 

Vs 

14.5 

2.03 


Parameter 

Value 

Standard Error 

C 

n 

o 

-0.000559 


c 

"6 

0.0506 

0.000596 

c 

n 

P 

-0.0769 

0.00289 

C 

n 

r 

-0.0878 

0.00716 

C 

"6 

a 

-0.000272 

0.000676 

C 

"6 

r 

-0.0505 

0.000553 
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Figure 2.- Configuration with droop outboard leading edge 
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Figure 5.- Model characteristics from stepwise regression technique. 
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Figure 6.- Longitudinal-force derivatives estimated from flight data 
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Figure 10. - Roll-rate derivatives estimated from flight data. 
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Figure 12,- Aileron derivatives estimated from flight data. 
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Figure 13.- Rudder derivatives estimated from flight data 
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Figure 18.- Comparison of sideslip derivatives determined from flight 
and wind tunnel- 
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Figure 25 •- Body system of axes. Arrows indicate positive direction 
of quantities. 
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